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Abstract 


This  paper  develops  methods  for  the  exact  computation  of  the 
distribution  of  the  maximum  flow  and  related  quantities  in  a 
planar  network  with  independent  and  exponentially  distributed  arc 
capacities.  A  continuous  time  Markov  chain  (CTMC)  with  upper 
triangular  rate  matrix  and  single  absorbing  state  is  constructed 
with  the  property  that  the  time  until  absorption  in  this 
absorbing  state  is  equal  to  the  value  of  the  maximum  flow  in  the 
network.  Recursive  algorithms  are  developed  for  computing  the 
distribution  and  moments  of  the  maximum  flow.  Algorithms  are 
also  presented  to  compute  the  probability  that  a  given  cut  is  the 
minimum  capacity  cut  in  the  network.  The  algorithms  are 
efficient  and  computationally  stable.  Distribution  of  the 
maximum  flow,  given  a  minimum  cut,  is  studied. 


Introduction 


1  . 

In  this  paper  we  investigate  directed  planar  stochastic 
networks  and  consider  the  problem  of  computing  the  probability 
distribution  of  the  maximum  flow  that  can  be  sent  from  the  source 
s  to  the  sink  t  without  violating  the  capacity  constraints  of  the 
arcs.  Among  several  situations  where  the  maximum  flow  problem  in 
stochastic  networks  is  encountered  in  practice,  we  describe  two. 

First,  consider  a  network  with  deterministic  arc  capacities. 
Assume  that  the  existing  flow  in  the  network  is  random,  thus  the 
excess  capacity  available  on  each  arc  is  random.  We  are 
interested  in  determining  how  much  additional  flow  can  be  sent 
from  the  source  to  the  sink  in  this  network.  Thus  we  are  faced 
with  the  maximum  flow  problem  in  a  stochastic  network  where  the 
arc  capacities  are  the  random  excess  capacities  in  the  original 
network.  Such  probabilistic  capacitated  networks  may  arise  in 
comm un i ca t i on  or  t r a n s po r t a t i on  applications  (see  Frank  and 
Frisch  [1971  ])  . 

Second,  consider  a  system  of  components  which  can  be 
represented  by  a  network.  Each  arc  of  the  network  represents  a 
component  in  the  system.  Each  component  can  be  either  "up"  or 
"down".  The  system  as  a  whole  is  "up"  if  and  only  if  there  is  a 
path  of  "up"  components  from  the  source  to  the  sink  in  the 
network  representation  of  the  system.  If  we  define  the  capacity 
of  an  arc  to  be  0  or  1  depending  upon  whether  the  component 
represented  by  the  arc  is  "down"  or  "up"  we  get  a  stochastic 
network.  The  probability  that  the  system  is  "up"  is  the 
probability  that  the  maximum  flow  in  the  stochastic  network  i3 
nonzero.  In  fact  the  actual  value  of  the  maximum  flow  can  be 
taken  to  indicate  the  "health"  of  the  system:  the  higher  the 
maximum  flow,  the  "healthier"  the  system. 


When  we  consider  multi-state  m ul t i- component  systems,  i.e. 


systems  where  the  components  as  well  as  the  system  may  exist  in 
more  than  two  states  of  degradation,  the  maximum  flow  concept 
provides  a  natural  extension  of  the  binary  case.  We  construct  a 
stochastic  network  where  the  capacity  of  an  arc  is  a  monotonic 
function  of  the  state  of  the  component  it  represents.  The  state 
of  the  system  is  then  defined  as  the  maximum  flow  from  the  source 
to  the  sink  in  this  network. 

The  value  of  the  maximum  flow  in  a  stochastic  network  is  a 
random  variable.  We  are  interested  in  its  distribution,  moments 
and  its  dependence  on  the  network  characteristics.  Frank  and 
Frisch  [1971]  have  offered  a  formulation  utilizing  transforms. 
The  procedures,  though  conceptually  simple,  are  impractical  from 
the  computational  point  of  view,  even  in  the  special  case  of 
normally  distributed  arc  capacities.  Evans  [1976]  has  addressed 
the  problem  for  a  network  with  integer->valued  ,  discretely 
distributed  arc  capacities. 

C  om  pr ehens  i  ve  procedures  have  been  developed  for  the  maximum 
flow  problem  when  arc  capacities  are  deterministic.  Ford  and 
Fulkerson  [  1  962  ]  developed  a  procedure  and  showed  that  the  value 
of  the  maximum  flow  from  the  source  s  to  the  sink  t  is  equal  to 
the  capacity  of  the  minimum  capacity  (s,t)  cut  in  the  network. 
When  the  network  is  (s,t)  planar,  i.e.,  it  has  a  planar  graph 
representati  on  such  that  an  arc  from  t  to  s  can  be  drawn  without 
violating  planarity,  one  can  draw  its  dual  network.  The  shortest 
path  problem  for  this  dual  network  can  then  be  solved  by  using 
one  of  the  standard  shortest  path  algorithms  (see  Dijkstra 
[1959],  Elmaghraby  [1977]).  The  duality  theory  asserts  that  the 


length  of  the  shortest  path  in  the  dual  network  is  equal  to  the 
capacity  of  the  minimum  capacity  cut  in  the  original  network,  and 
thus  provides  the  value  of  the  maximum  flow.  Alternate 
algorithms  for  the  maximum  flow  in  (s,t)  planar  networks,  which 
avoid  constructing  the  dual,  are  described  in  Itai  and  Shiloach 
[1979]. 

In  this  paper,  we  consider  the  maximum  flow  problem  in  (3,t) 
planar  networks  where  arc  capacities  are  independent 
exponentially  distributed  random  variables.  When  the  capacities 
are  not  exponentially  distributed  one  can  approximate  their 
actual  distributions  by  phase  type  distributions.  The  methods 
developed  in  this  paper  can  be  directly  extended  to  the  case  when 
the  arc  capacities  have  phase  type  distributions.  We  do  not 
discuss  this  extension  in  this  paper  since  the  most  significant 
contribution  of  this  paper  is  the  construction  of  the  stochastic 
processes  which  facilitate  the  solution  of  the  maximum  flow 
problem . 

In  Section  2  we  establish  the  network  terminology  related 
to  (s,t)  planar  networks.  In  particular  we  introduce  the 
"topmost  path  first"  order  on  the  set  of  all  (s,t)  paths  in  the 
network.  In  Section  3  we  present  an  algorithm  for  the  maximum 
flow  problem  in  (s,t)  planar  networks  with  deterministic  arc 
capacities,  which  is  subsequently  used  in  the  development  of  our 
al gor i thms . 

In  Section  4  we  construct  a  finite  state  space  continuous 
time  Harkov  chain  (CTMC )  with  the  following  property:  the  time 
until  absorption  in  the  absorbing  state  in  this  Markov  chain  is 
equal  to  the  value  of  the  maximum  flow  in  the  network.  The  state 


apace  is  shown  to  be  the  set  of  (s,t)  paths  in  the  network 
augmented  by  an  absorbing  state.  When  these  paths  are  ordered 
according  to  the  "topmost  path  first"  order,  the  rate  matrix  of 
the  CTMC  becomes  upper  triangular.  This  structure  yields 
computationally  simple  algorithms  to  compute  the  distribution  and 
moments  of  the  maximum  flow.  These  algorithms  are  described  in 
Section  5. 

In  Section  6  we  introduce  the  criticality  indicator  of  a 
minimal  (s,t)  cut,  which  takes  the  value  1  if  the  cut  is  a 
minimum  capacity  cut  in  the  network,  and  0  otherwise.  The 
criticality  index  of  a  cut  is  defined  to  be  the  probability  that 
the  criticality  indicator  of  the  cut  is  1.  We  study  the  joint 
distribution  of  the  criticality  indicator  of  a  cut  and  the 
maximum  flow  and  develop  algorithms  to  compute  the  criticality 
index  of  a  cut. 

In  Section  7  we  document  the  computational  results  using 
four  networks  as  examples.  The  conclusions  are  presented  in 


Section  8. 


2. 


Planar  Networks 


Let  G  =  ( V , A )  be  a  network  with  vertex  set  V  ,  arc  set  A, 
source  node  s,  sink  node  t,  and  N  paths.  It  is  assumed  that  G  is 
(s,t)  planar  and  that  one  planar  representation  of  G  is  fixed  in 
advance . 

An  arc  eeA  is  incident  on  a  node  u  if  u  is  the  starting  or 
ending  node  of  e.  Let  E(u)  denote  the  set  of  arcs  incident  on  u. 
Each  arc  eeE(u)  induces  a  partial  order,  on  E(u)  as  follows.  Arc 
e^  is  said  to  be  bef  ore  e2  i_n  the  partial  order  i  nduced  by  e  on 
E(u)  if  ei,  e2eE(u)  and  ei  appears  before  e2  in  a  clockwise  sweep 
around  u  starting  from  e. 

A  path  in  G  is  described  by  the  sequence  of  nodes  it  visits. 
A  path  is  called  simple  if  it  does  not  visit  any  node  more  than 
once.  Let  Pi  *  ( v 0 - s  ,  v i ,  .  .  . , vm - t  )  and  P2”(uo”3  »  U1  »  •  •  •  * un«t )  be 
two  simple  (s,t)  paths  in  G.  Assume  that  an  imaginary  arc  (t,s) 
is  drawn  on  the  planar  representation  of  G  without  violating  the 
planarity.  Let  v*i«t.  The  path  Pi  lies  above  P2  (or, 
equivalently,  P2  lies  below  Pi)  if  there  exists  an  index  r  such 
that  vi»ui  for  i«1,2,...,r  and  vr  +  i  ^  ur  +  i  and  the  arc  (vr,vr+i> 
is  before  (ur,ur  +  i)  in  the  partial  order  induced  by  (vr^i  ,vr)  on 
E  ( vr  )  . 

The  relation  "lies  above"  defines  a  partial  order  on  the  set 
of  all  (s,t)  path3  in  G.  Let  L  be  the  ordered  set  of  paths 
{ Pi  , P2 »  •  •  • » ?N 1  3uch  that  if  Pj  lies  above  Pj  then  i  <  j.  The 
paths  in  L  are  said  to  be  in  the  topmost  path  first  (TPF)  order. 
An  algorithm  to  list  the  paths  in  the  TPF  order  is  described 
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below.  It  1 s  a  modified  version  of  the  backtracking  algorithm 
described  in  Reed  and  Tarjan  [1975]. 


Step  1 
Step  2 


Algorithm  to  List  Paths  i n  TPF  Order : 

Step  0.  Set  v^i«t. 

All  nodes  unlabelled. 

All  arcs  unscanned. 

Set  i=0,  vg-s.  Label  s. 

Let  e  be  the  first  unscanned  arc  leaving  vA  that  appears 
in  the  clockwise  sweep  around  v^  starting  from  (v^^  ,Vj  )  such 
that  head  node  of  e  is  unlabelled. 

If  no  such  e  exists,  go  to  Step  4 . 

Step  3:  Set  vi+i  »  head  node  of  e,  label  vi+1.  Scan  e. 

If  v<  »  t,  go  to  Step  5. 

else,  set  i=l+1  and  go  to  Step  2. 

Step  .  If  i»0,  stop.  (All  paths  are  listed). 

Else,  unscan  all  arcs  leaving  v^,  unlabel  v^. 

Set  i=i-1  and  go  to  Step  2. 

Step  5.  An  (s,t)  path  of  length  i+1  is  found. 

List  the  nodes  vq,vi . v^. 

Unlabel  t,  set  i-i-1  and  go  to  Step  2.  □ 


Choice  of  e  in  Step  2  guarantees  that  the  paths  are  generated  in 
the  TPF  order. 

Let  P i - ( v g-s  ,  v i , .  .  . , vm« t )  and  P2 “ ( uo"3 »ui . un« t )  be  two 

simple  (s,t)  paths  in  G.  The  path  P^  lies  completely  above  P2 
(or,  equivalently,  P2  lies  completely  below  Pi )  if  vj-uj  and  vi+i 
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f  Uj+1  implies  that  (vitvi+1)  is  before  (uj,uj  +  1)  in  the  partial 
order  induced  by  (vjHi»Vi>  on  E(vi).  For  each  P  in  L  and  an  arc 
e  on  P,  define  an  alternate  path  P(e)  to  be  the  first  path  in  L 
that  does  not  use  arc  e  and  lies  completely  below  P.  If  no  such 
path  exists  we  say  P(e)  -  $.  An  algorithm  to  construct  P(e) 
given  P  and  e  is  described  in  Itai  and  Shiloach  [1979]. 

We  end  this  section  with  an  example  which  illustrates  the 
concepts  defined  here.  Consider  the  planar  network  of  six  nodes 
and  nine  arcs  displayed  in  Figure  1 .  Node  1  is  the  source  and 
node  5  is  the  sink.  Table  I  lists  the  paths  in  this  network 
in  TPF  order.  Path  ( 1 ,3*^*5)  lies  below  path  (1,2, 3,6, 5)  but 
does  not  lie  completely  below.  Path  (1  ,3,4,5)  lies  completely 
below  path  (1,2, 3,4,5).  The  alternate  path  for  path  (1,2, 4, 6, 5) 
is  (  1  ,  2  ,  3  ,  4 , 6 , 5  )  if  arc  e  is  (2,4).  There  is  no  alternate  path 
for  path  (1,2, 4, 6, 5)  if  arc  e  is  (6,5). 


3.  Maximum  Flow  i n  Deterministic  Planar  Networks 

In  this  section  we  describe  an  algorithm  to  find  the  maximum 
flow  from  s  to  t  in  a  planar  network.  The  algorithm  is  suggested 
in  Ford  and  Fulkerson  [1956],  developed  in  Berge  and 
G ho u i 1  a- H o ur i  [196  2  ]  and  its  time  complexity  is  reduced  to 
0(  |v|  log  |v|  )  by  Itai  and  Shiloach  [  1  97  9  ].  The  algorithm  is 
described  below  since  it  is  the  starting  point  for  our  analysis 
in  the  stochastic  case. 

Let  c(e)  be  the  non  negative  deterministic  flow  capacity  of 

arc  e.  Let  f*  be  the  maximum  flow  that  can  go  from  s  to  t.  As 

before  let  L  be  the  list  of  (s,t)  paths  in  the  TPF  order.  Define 

the  capacity  of  a  path  PeL  as 

c( P)  =  min  {c  (e )  1  .  (31) 

e  o  n  P 

Maximum  Flow  Algorithm: 

Step  1 .  Let  P  be  the  first  path  in  L.  Set  f*  =  0. 

Step  2.  f*  =  f*  +  c ( P ) 

c(e)  =  c(e)  -  c ( P )  for  all  eeP. 

Step  3 .  Let  P*  be  the  next  path  in  L  such  that  c(P')  >  0. 

If  no  such  path  exists,  go  to  Step 

Else,  set  P  <-  P  ’  and  go  to  Step  2. 

Step  .  Stop,  f*  is  the  maximum  flow.  □ 

Remark  1  .  This  algorithm  is  called  the  "path  filling"  algorithm, 
since  it  starts  with  the  first  path  in  L  and  keeps  filling  them 
to  capacity  until  no  more  path  can  be  found.  The  proof  that  this 
algorithm  produces  the  maximum  flow  is  given  by  both  Berge  and 
Ghouila  Houri  [1962]  and  Itai  and  Shiloach  [1979]. 


Remark  2  .  In  Step  2,  let  e  be  that  arc  on  P  for  which  c(P)  = 
c(e).  If  this  e  is  unique,  then  it  can  be  shown  that  P'  is  the 
alternate  path  P(e)  as  defined  in  Section  2. 

Remark  3.  The  path  filling  algorithm  works  if  the  paths  are  used 
in  TPF  (or  the  reverse  of  TPF)  order.  An  arbitrary  order  may  not 
work,  i.e.  all  the  paths  may  get  saturated  while  the  flow  is  not 
maximum . 

Remark  4  .  The  TPF  order  is  not  defined  for  a  non  planar  network 
In  such  a  case  it  is  possible  to  construct  an  example  where  the 
path  filling  algorithm  fails  no  matter  what  ordering  of  paths  is 


-10- 


*1 .  Maximum  Flow  in  Stochastic  Planar  Networks 


Let  G  »  ( V ,  A )  be  the  (s,t)  planar  network  with  source  s  and 
sink  t  and  L  be  the  ordered  set  of  (s,t)  paths  in  G  in  TPF  order. 
Let  C(e)  be  a  non-negative  random  variable  representing  the 
capacity  of  arc  e.  Let  T  be  the  maximum  flow  that  can  be  sent 
from  s  to  t.  We  are  interested  in  the  distribution  and  moments 
of  T.  To  analyze  this  problem  we  construct  a  stochastic  process 
with  a  single  absorbing  state  so  that  the  time  until  absorption 
in  this  state  is  equal  to  the  maximum  flow,  T,  from  s  to  t. 

Visualize  the  network  G  as  a  network  of  pipelines,  with  arc 
e  capable  of  handling  C(e)  gallons  of  flow  per  unit  time. 
Suppose  that  there  i3  a  pump  at  node  s  that  forces  fluid  into  the 
network  at  the  rate  of  t  gallons  per  unit  time  at  time  t.  Thus 
the  flow  of  the  pump  increases  linearly  with  time.  The  nodes  act 
as  multiple-valves  controlled  by  a  complex  mechanism  which 
operates  as  follows:  At  time  t  »  0 ,  the  first  path  in  L,  call  it 
P,  is  open  for  the  fluid  flow;  all  other  paths  are  blocked.  When 
this  path  fills  up  due  to  capacity  restriction  of  arc  e  on  P,  the 
flow  is  then  diverted  to  the  alternate  path  P(e).  This  process 
continues  as  long  as  there  is  an  unfilled  alternate  path 
available.  When  all  paths  are  filled  the  pump  flow  cannot  be 
increased  any  further. 

Define  X(t)  to  be  the  path  in  which  flow  is  increasing  at 
time  t.  If  at  time  t  the  network  is  saturated,  i.e.  the  flow  is 
not  increasing,  X(t)  is  defined  to  be  $.  X(t)  represents  the 
state  of  the  system  at  time  t.  Note  that  the  network  handles  a 


flow  of  t  gallons/unit  time  at  time  t.  From  the  construction  of 


X ( t )  and  the  path  filling  algorithm  of  Section  3,  it  is  clear 


that 

T  -  min  {x(t )  -  ♦} .  (*».  1 ) 

t20 

From  now  on  we  assume  the  following: 

A 1  .  {C(e),  eeA}  is  a  set  of  independent  random  variables. 

A2.  C(e)  is  exponentially  distributed  with  parameter 
u(e)>0  (mean  -  1/p(e)). 

Theorem  1  .  Under  assumptions  A1  and  A2,  (X(t),  1 2;  0 }  is  a  CTMC 

with  state  space  L*  -  Lu  }  $  }  and  rate  matrix  Q  *  [q(P,P')], 
P  ,  P  '  eL*  ,  where 

q(P.P')  -  I  u(e)  if  ?'{  P;  P.P’eL 

e  on  P:  P(e)  «  P  ' 

-  -I  p(e)  if  P  -  P'eL  (4.2) 

e  on  P 

»  0  i  f  P  =  <j> . 

Proof:  The  fact  that  {X(t),t20}  is  a  CTMC  follows  from 

i 

assumptions  A1  and  A2.  Now  suppose  that  X(t)  -  PeL.  If  the  arc 
e  on  P  saturates,  which  happens  at  rate  y(e),  the  state  changes 
to  P(e).  If  P(e)  does  not  exist,  the  state  changes  to  <|> .  If 
X  ( t )  «  <j>  it  remains  41  from  then  on.  The  equations  for  q(P,P’) 
follow  from  these  observations.  a 

Remark  5  .  From  the  construction  of  fx(t),t*o}  it  is  obvious  that 
the  process  does  not  visit  the  same  state  twice,  since  once  a 
path  is  saturated  it  stays  saturated.  Hence  the  rate  matrix  is 
upper  triangular.  Thi3  structural  property  is  critical  in 


developing  simple  algorithms  for  computing  the  distribution  and 
moments  of  T. 

Example.  Consider  the  example  network  of  Figure  1.  The  state 
space  for  the  CTMC  representation  for  this  Markov  chain  is  partly 
given  in  Table  I.  (Table  I  gives  L,  the  state  space  is  L*  - 
Lul<J>}).  The  rate  matrix  for  this  problem  is  given  in  Table  II. 


5.  Distribution  and  Moments  of  Maximum  Flow 


In  this  section  we  develop  algorithms  to  compute  the 
distribution  and  moments  of  the  maximum  flow.  It  is  possible  to 
develop  two  types  of  algorithms,  called  backward  and  forward 
algorithms.  We  present  both  versions,  but  we  illustrate  only  the 
backward  algorithms  by  means  of  an  example.  Also,  from  now  on  we 
identify  the  set  L  with  the  set  {l,2,...,N},  i.e.  X(t)  -  i  means 
that  the  state  of  the  system  at  time  t  is  given  by  the  i-th  path 
in  L.  The  state  <j>  will  be  identified  by  the  number  N  +  1.  Thus 
equation  (4.1)  becomes 

T  -  min  { t SO :  X(t)  -  N+1  |  X(0)  -  if.  (5.1) 

We  are  interested  in  the  distri bution  of  maximum  flow 

F(t)  -  P ( TSt )  (5.2) 

and  the  k-th  moment 

m(k)  -  E(Tk)  .  (5.3) 


5 . 1  Distribution  of  Maximum  Flow 
Backward  Algor  1 thm :  Define 

pt(t)  -  P|x(t)  -  N  + 1  |  X(0)  »  i  }  ,  ISiSN  +  1.  (5.4) 

Then  F(t)  «  pi(t).  The  differential  equations  for  p^(t)  are 


j,  0  otherwise.  Due  to  the  upper 


where  6i,j  *  1  if  i  - 

triangular  nature  of  Q,  equations  (5.5)  can  be  solved  in  a 
backward  manner,  starting  with  pN+j(t)  -  1  for  t  £  0,  and 

computing  pjj  ( t ) . p-|(t)  -  F(t)  in  that  order. 

Example .  For  network  1,  equations  (5.5)  become 

Pg(t)  -  0, 

Pg(t)  -  -(u2  +  y^  +  y^)Pg(t  )  +  ( u2  +  Ug  +  Ug)Pg(t  )  , 

pj(  t)  -  “  (  m2  +  g  +  M  g  )  P^  (  t  )  +  (U5  +  Wg)Pg(t)  +  (  U2  +  >J9^P9^  t  ^  ' 

:  (5.6) 

P^(t)  -  -(m1  +ul4  +  u5  +  u7)p^(t )  +  UyP^  (t )  +  \i  ^Pg(t)  +  (  u1 +y  )  Pg  ( t ) , 

P  2  (  t  )  »  -  (  M1 +U  g  +  U  g  +  Mg)  P2  (  1  ^  +  ( Wg+Ug)P5(t )  +  ^^7^  +  PqPg^t)* 

pj(t)  -  -'(u1+u3  +  m7)p1  (t)  +u7P2(t)  ♦  Ji^U)  +  u,P6(t)- 

Forward  Algor i thm :  Define 

p.  (t)  -  P { X ( t  )  -  j|X(0)  -  1},  1  S  j  S  N  ♦  1.  (5.7) 

J 

Then  F(t)  -pfj  +  i(t).  The  differential  equations  for  Pj(t) 
are  given  by 


Pj(t) 


I  P^tXi 
iSj  J 


1  $  j  S  N  ♦  1 


(5.8) 


p.(0)  -  6jt1. 

Again  since  Q  is  upper  triangular,  the  above  equations  can  be 
solved  in  a  forward  manner,  starting  with  p^(t)  -  exp(q^t),  t  2  0, 
and  computing  P2 ( t  ),...,  pN + j ( t  )  »  F(t)  in  that  order. 

Numerical  Evaluation  of  F ( t  ) :  U n i f ormi zat i on  technique  described 
in  Ross  [1983]  is  used  to  numerically  evaluate  F(t).  The  details 


of  this  technique  are  given  in  Kulkarni  and  Adlakha  [1984].  Here 


we  give  only  the  final  result.  Define 

q  -  max  j-^q.  }  (5.9) 

1SiSN+1 

and  let 

q* j  -  6ij  +  q^/q,  1  S  1.  J  S  N  +  1  (5.10) 

#  # 

Let  Q  -  [q  ]  and  let  an  be  the  (1,N+1)th  element  of  the  nth 

* 

power  of  Q  .  It  is  easy  to  see  that  a  's  increase  to  1  as  n  +  ®. 

n 

Then  we  have 

CD 

F(t)  “  l  «  e  Qt(qt)n/n!.  (5.11) 

n-0 

In  numerical  calculations  we  choose  an  e  >  0  and  an  M  such  that 
2  1«e.  Then  compute 


FM(t) 


M  _ 

l  (1-a  )e  Qt(qt)n/n! 
n  =  0 


(5.12) 


and 


F„<t) 


M 

l 

n  =  0 


( a* 


'an)e 


qt(qt)n/n! 


(5.13) 


It  can  be  shown  easily  that  F^(t)  and  ^^(t)  provide  lower  and  upper 
bounds  for  F(t)  within  e.  It  is  a  straightforward  matter  to 
compute  an's  sequentially.  Thus  uniformization  is  a 
computationally  stable  and  rapid  method  of  computing  F(t). 
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5 . 2  Moments  of  Maximum  Flow 

We  now  describe  the  backward  and  forward  algorithms  to 
compute  m(k),  the  k->th  moment  of  T. 

Backward  Algorithm;  Define 

Ti  -  min  (t*0:  X(t)  -  N  +  1  |  X(0)  -  i|,  1  S  i  S  N  +  1  (5.1*1) 


and 

m t ( k )  -  E(T*)  ,  1  S  i  S  N  ♦  1 .  (5.15) 

As  { X ( t  )  ,  t20}  gets  absorbed  in  the  state  N  +  1  with  probability 
1,  nijjo)  -  1  for  TSiSN  +  1.  Also,  for  k  2  1,  mN+i(k)  -  0. 
It  can  be  shown  that 

m  (k)  -  [  km  ( k->  1  )  +  £  q..m  (k)]/q  (5.16) 

J>i  J  3 


where  -  n  q  ±  A  .  We  have  m(k)  -  m^k).  To  compute  m(k)  one 
needs  to  compute  m1(r)  for  r  -  1,2,...,k;  i  -  N  ♦  1,  N,...,2,1  in 
that  order. 

Example .  For  network  1  equations  (5.16)  become 
m  q  (  k  )  «  0  , 


mg(k  ) 
m7  ( k  ) 


[  kmg ( k- 1 ) 
[km7 (k- 1 ) 


♦  (  u2  +  u6  +  Ug)mg(k ) ]/ (u2  +  y6  +  y9)  , 

+  (u5  +  ug)m8(k)  +  (  u2  +  ug)m9(k  )  ]/ (n2  +  y5  +  ,j8  +  n9)  , 

(5.17) 


(k  ) 
m2(k  ) 
i"1  (k  ) 


[km j (k-  1 ) 
[ km2  (k-  1  ) 
[km1 (k-  1  ) 


+  u7m2j(k)  ♦  n5m5(k)  ♦  (  y,  +  y  4  )mg  (  k  )  ]/ (  +  y  ^  +  yg  +  y  ) 

+  (u3  +  u8)m5(k )  +  y  1  m 7 ( k )  ♦  y 9mg ( k )  ]/ ( y }  * y ^ ♦ v g  +  y g ) 

+  Vjm2 ( k )  ♦  y^U)  ♦  y 1 mg ( k  )  ]/ ( M 1 ♦ u  ♦ M ^ ) . 


Forward  Algorithm:  Define 


inf jt*0:X(t)  -  J  |  X(0)  -  1  } 


(5.18) 


where  the  infimum  over  an  empty  set  is  defined  to  be  +  *.  Let 


mj (k )  -  E( Tj ) , 


(5.19) 


It  is  possible  to  have  m j ( 0 )  ■  P(Tj  <  »)  <  1  for  1  £  J  S  N. 
It  can  be  shown  that 


^  ^  L/  m  V* 

m  (k)  -  l  [  l  (k! /r!  )  (m.  (r  )/ q  )](q11/Qi) 
J  i < J  r-0  J 


(5.20) 


ral(k)  *  5k,0 


k*0,  1  S J  SN ♦ 1  . 


These  equations  are  solved  recursively  in  a  forward  manner.  It 
is  clear  that  equations  (5.16)  can  be  solved  more  efficiently 


than  equations  (5.20). 


6.  Joint  Performance  and  Criticality  Index 


Let  Yc A  be  a  prespecified  (s,t)  cut.  Assume  that  Y  is  a 
minimal  cut,  i.e.,  no  proper  subset  of  Y  is  a  (s,t)  cut.  Then  it 
is  known  that  there  is  a  set  BcV  such  that 

Y  -  C(B,B)  -  { ( u  ,  v  )  e A  :  ueB.vel}.  (6.1) 

Define  the  capacity  of  a  cut  to  be  the  sum  of  the  capacities  of 
its  constituent  arcs.  Let 

I ( Y )  -  1  if  Y  is  minimum  capacity  (s,t)  cut  in  G 

=  0  otherwise.  (6.2) 

I  ( Y )  is  called  the  criticality  indicator  of  the  cut  Y.  Define 
the  criticality  index  of  the  cut  Y  as 

R ( Y  )  -  P ( I ( Y  )  «  1  )  .  (6.3) 

In  this  section  we  study  the  joint  distribution  of  I(Y)  and  the 
maximum  flow  as  well  as  the  criticality  index  R ( Y ) .  The  key  to 
the  development  of  these  algorithms  is  the  following 
characterization  of  minimum  capacity  cut  from  Ford  and  Fulkerson 
[1  962]: 

A  cut  Y  «  C  ( B ,  B )  i_s  a  minimum  capac  i  ty  (s,t)  cut  if  and 
only  i f  the  arcs  in  Y  are  saturated  while  the  arcs  i n  C(B,B) 
are  empty  when  the  flow  from  s  to;  t  i_s  maximum. 

A  (s,t)  path  in  G  is  called  Y-admlsslble  if  it  intersects  Y 
in  one  and  only  one  arc.  It  is  clear  that  a  path  is  Y-admissible 
if  and  only  if  it  has  no  arcs  from  C(B,B).  It  can  be  easily 
seen  that  the  first  and  the  last  paths  in  L  are  Y-admissible  for 
any  minimal  cut  Y.  Let  L  ( Y  )  be  the  ordered  set  of  Y-admissible 
paths  obtained  by  deleting  all  non  Y-admissible  paths  from  L. 
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We  now  modify  the  pipe  filling  process  described  in  Section 
4  as  follows:  Consider  the  following  cases  when  the  current 
path,  say  P,  fills  up  due  to  the  capacity  restriction  of  arc  e  on 
P: 


(1 ) 

eeY  , 

P  ( e ) 

» 

♦ » 

(2) 

e£Y, 

P(e) 

a 

$  > 

(3) 

eeY  , 

P  ( e ) 

* 

$ » 

P(e) 

i  s 

Y-admissi ble , 

(4) 

e£Y , 

P(e ) 

* 

+  . 

P(e  ) 

i  s 

Y- admi ssi bl e  , 

YnP 

-  YnP(e), 

(5) 

e£Y  , 

P(e) 

# 

♦  • 

P(e) 

i  s 

Y-admissible, 

YnP 

*  Yn  P  ( e  )  , 

(6) 

e£Y  , 

P  ( e ) 

# 

4>  • 

P  ( e ) 

i  s 

not  Y-admissible. 

It  can  be  shown  that  the  case  "eeY,  P(e)  *  $ ,  P(e)  is  not 
Y- adm i ss i bl e"  cannot  occur  if  Y  is  a  minimal  cut.  The  original 
pipe  filling  process  described  in  Section  4  terminates  if  cases  1 
or  2  are  encountered.  The  modified  process  terminates  if  cases 
1,  2,  5  or  6  are  encountered.  If  the  modified  process  terminates 
due  to  case  1,  It  will  be  called  a  "good"  termination;  else  it 
will  be  called  a  "bad"  termination.  in  cases  3  or  4  the  process 
continues  by  diverting  the  flow  to  the  path  P(e). 

Let  Xy ( t )  represent  the  path, in  which  flow  is  increasing  at 
time  t  in  the  modified  process.  if  the  flow  is  not  increasing  in 
any  path  at  time  t,  the  process  must  have  terminated  before  t. 
If  the  process  ends  with  a  "good"  ("bad")  termination  at  time  t, 
Xy(t)  is  defined  to  be  <j>(4>')  for  t  2  t.  It  is  clear  that  as  long 
as  X  y  ( t  )  e  L  ,  it  is  in  fact  in  L  ( Y  )  ,  since  Xy(0)  is  a  Y-admissible 
path.  Thus  {Xy(t),t20}  is  a  stochastic  process  with  state  space 

L»(Y)  -  L(Y)uU,*'  }.  (6.4) 


The  following  theorem  brings  out  the  significance  of  the  process 
i X Y (t ) ,t*0} . 


Theorem  2 .  Let  {Xy(t),t20}  be  the  modified  pipe  filling  process 
on  L»(Y)  with  XY(0)  -  first  path  in  L ( Y )  .  Let 

T ( Y )  -  min  {tSO:  X  (t)  -  *}.  (6.5) 

Then 

P(T(Y)St)  -  P ( I ( Y )  -  1  and  capacity  of  YSt).  (6.6) 

Proof :  We  shall  show  that 

( X  y  ( t )  -  $ }  »  1 1  ( Y  )  =  1  and  capacity  of  Y  S  t}.  (6.7) 

(->)  Suppose  X  y  ( t  )  -  4> .  Then,  from  the  construction  of  the 
modified  pipe  filling  process,  maximum  flow  of  value  S  t  has  been 
achieved;  no  arcs  of  C(B,B)  have  been  used  and  the  arcs  in  Y  are 
full.  Then  from  the  characterization  of  max-flow  by  min-cut,  the 
cut  Y  is  the  minimum  capacity  cut  and  its  capacity  (which  is  same 
as  the  maximum  flow)  is  S  t.  Thus  I(Y)  »  1  and  capacity  of  Y  S  t . 
(<-)  Suppose  I  ( Y  )  -  1  and  capacity  of  Y  S  t.  Then  from  the 
characterization  o'1  max-flow  by  min-cut  the  arcs  in  C(B,B)  are 
empty,  the  arcs  in  Y  are  saturated  and  the  maximum  flow  <  t. 
Hence,  XY(t)  must  have  stayed  in  L(Y)  until  it  hit  <p  before  t. 
Thus  X Y  ( t )  =  <p  . 

As  <j>  is  an  absorbing  state 

P ( T ( Y  )  St)  -  P(XY(t)  -  $) .  (6.8) 

The  result  now  follows  from  equations  (6.7)  and  (6.8).  n 

A3  <(>  '  is  also  an  absorbing  state,  T(Y)  is  an  incomplete 
random  variable.  In  fact 

R(Y)  =  P ( I ( Y  )  -  1)  -  P ( T ( Y  )  <  »).  (6.9) 

Thus,  the  criticality  index  of  a  cut  Y  is  given  by  the 
probability  that  the  process  |xY(t),  1 2  0 }  gets  absorbed  in  <j>. 


Note  that  Theorem  2  is  true  regardless  of  the  assumptions  A1  and 
A2.  If  A 1  and  A2  hold,  }XY(t),  t£0}  oan  be  shown  to  be  a  CTMC  as 


proved  in  the  following  theorem.  First,  we  introduce  some 
notation. 


If  X y ( t )  «  PeL(Y)  and  arc  e  on  P  saturates,  a  transition 

takes  place  to  some  other  path  P'eL(Y),  41  or  4 '  .  This  transition 
is  denoted  by  (P,e)  and  belongs  to  one  of  the  six  cases  described 
earlier.  For  PeL(Y),  define 

M ( P , P ' )  «  {e  on  P:  P'  =  P(e),  the  transition 

(P,e)  belongs  to  case  3  or  4},  (6.10) 

M(P,4>)  =>  {e  on  P:  the  transition  (P,e) 

belongs  to  case  l},  (6.11) 

M  (  P ,  (j> '  )  -  {e  on  P:  the  transition  (P,e) 

belongs  to  case  2,  5  or  6).  (6.12) 

Theorem  3.  Under  assumptions  A1  and  A2,  jXy(t),  t£o}  is  a  CTMC 
on  L*(Y)  with  the  rate  matrix  Q(Y)  *  [qy^.P')]  (  P  ,  P  '  e  L  *  (  Y  )  )  , 

defined  as  follows 


q  y  (  p  »  p  ' )  - 

l  u(e ) 

eeM ( P  ,  P  '  ) 

if 

P,P’ 

eL ( Y  )  ,  P*P’ 

s 

l  u(e ) 

eeM(P,<j>) 

if 

P  '  = 

4>,  PeL(Y) 

s 

l  u(e) 

eeM ( P , 4  '  ) 

if 

P  '  - 

<t>',  PeL(Y) 

- 

~l  u(e) 

e  on  P 

if 

P '  ■= 

P  eL  ( Y  ) 

- 

0 

if 

P  = 

4  or  4  '  . 

Proof:  Follows  from  the 

construction  of 

{  X  Y ( t )  ,  t>0  } 

def i ni ti ons 

of  M ( P , P ' )  ,  M  (  P 

,  4> )  and  M  (  P  , 

<p'  )  . 

(6.13) 


and  the 
□ 


Let  the  paths  in  L(Y)  be  numbered  1,2 . m  and  the  states  <j> 

and  be  identified  by  the  numbers  m+1  and  m  +  2  respectively. 


Then  from  Theorem  2  we  get 

P ( I ( Y  )  -  1  and  capacity  of  Y  S  t) 

-P(I(Y)  -  1  and  maximum  flow  £  t) 

»P(Xy(t)  -  m+1  |  Xy(0)  »  1).  ( 6 . 1 4 ) 

Thus  the  joint  distribution  of  the  maximum  flow  and  the 
criticality  indicator  of  the  cut  Y  can  be  obtained  as  the 
distribution  of  time  until  absorption  in  the  state  m+1  of  the 
modified  CTMC  {Xy(t),t20}  starting  from  state  1.  The 
unif ormi zation  technique  of  Section  5  can  be  used  to  compute  this 
distribution  numerically. 

The  criticality  index  of  the  cut  Y  is  given  by 
R(Y)  =  P ( X y ( t  )  =  m  +  1  for  some  t  £  0  |  Xy(0)  =  1).  (6.15) 

Thus  the  criticality  index  can  be  computed  as  the  probability 
that  the  process  1  X  Y  (  t  )  ,  t  £  o},  starting  in  state  1,  gets 
absorbed  in  state  <p  .  The  backward  and  forward  algorithms  for 
doing  this  are  described  below.  We  write  qjj(Y)  for  qy ( P , P ' ) 
where  P  and  P'  are  identified  by  numbers  i  and  j  respectively. 
Backward  A1 gor i thm :  Define 

vi  “  P { X y ( t  )  =  m+1  for  some  t20  |  Xy(0)  -  i  }  ,  1SiSm+1.  (6.16) 

Then 

v  -  l  v  q  (Y)/q.(Y)  ISiSm  (6.17) 

j  >  i  J  J 

where  q^(Y)  -  -  q  i  i  ( Y  )  and  vm  +  -j  =  1,  vm  +  2  -  0.  These  quantities 
can  be  computed  in  a  backward  manner  to  yield  R(Y)  ■  v  ■)  . 

Example .  Consider  network  1  and  let  Y  be  the  minimal  cut 
[3.5,91.  For  this  case,  B  -  {l .2,3,6}  and  C(B,B)  =  { 8 } .  The 


set  of  Y  a  dm  i  s  s  i  bl  e  paths,  L(Y),  is  shown  in  Table  III.  The  Q(Y) 
matrix  is  shown  in  Table  IV.  Equations  (6.17)  for  this  problem 
are  as  follows: 


0, 


v  6  =  1* 


U9V6/  (,J2  +  V16',’U9)  ’ 

(u2  +  u5  +  p? )  , 


(6.18) 


v3  “  ((v1! +V‘j,)v5  +  1j9v6)/ (u1 +W4  +  U6  +  P9)  » 

V2  =  (w5V3  +  (vi1  +m4)vii)/ (ui  +Ui4  +  M5  +  M7)  ’ 

V1  *  m2v2/(p1+u3+u7). 

Solving  these  in  a  backward  fashion  we  get  =  p|the  cut  {3 
is  the  minimum  capacity  cut}. 

Forwar d  Algorithm :  Define 

Vj  =  p{xY(t)  =  j  for  some  tSO  |  Xy(0)-l},  lSJSm+1.  (6.19) 


Then  R ( Y )  =  vm+i*  The  equations  for  are  given  by 


vi  -  l  (q.  .  (Y  )/q  (Y)  )v  , 
J  i<J 


(6.20) 


which  can  be  solved  in  a  forward  manner  starting  with  v^  =  1 
computing  v2 . vm+i  *  R(^)  in  that  order. 


,5,9} 


and 


7 .  Computational  Results . 

In  this  section  we  document  the  computational  results  using 
the  networks  shown  in  Figures  1,2,3  and  H.  The  networks  are 
assumed  to  satisfy  assumptions  A1  and  A2  of  Section  4.  The 
backward  algorithms  developed  in  Sections  5  and  6  are  implemented 
in  Fortran  77. 

Table  V  gives  the  following  network  size  descriptors  for  the 
four  networks : 

(1)  Number  of  nodes, 

( 2)  Number  of  arcs  , 

(3)  Number  of  (s,t)  paths, 

( H)  Number  of  (s,t)  cuts,  and 

(5)  Value  of  the  un i f or m i za t i on  constant  q  (defined  by 
equation  (5.9)). 

Table  VI  gives  the  mean  (p)  and  standard  deviation  (o)  of 
the  maximum  flow  in  the  four  networks.  The  distribution  of  the 
maximum  flow  is  computed  using  the  uniformization  technique 
described  in  Section  5.  Thus,  an  integer  M  was  chosen  so  that  004 
>  1-10"5.  The  values  of  M  and  are  also  given  in  Table  VI. 

Note  that  one  does  not  necessarily  need  larger  M  for  larger 
networks . 

The  cumulative  distribution  function  F(t)  is  approximated  by 
FM(t)  of  equation  (5.12).  The  distribution  of  the  normalized 
maximum  flow,  i.e.  (maximum  flow-y)/o,  is  tabulated  in  Table  VII. 
An  interesting  feature  of  this  table  is  that  the  distribution  of 
the  normalized  maximum  flow  seems  to  be  more  or  less  the  same  for 
all  the  four  networks.  No  theoretical  insight  into  this 


phenomenon  is  currently  available,  other  than  the  intuitive 


feeling  that  for  "large"  networks  the  distribution  should  be 
close  to  normal  . 

The  conditional  performance  of  the  four  networks  is 
tabulated  in  the  next  four  tables.  For  each  network  the  tables 
list  the  following: 

(1)  R ( Y  )  ,  the  probability  that  Y  is  the  minimum  capacity 

cut , 

(2)  p ( Y  )  ,  the  mean  of  the  maximum  flow  given  that  Y  is  the 
minimum  capacity  cut, 

(3)  o(Y),  the  standard  deviation  of  the  maximum  flow  given 
that  Y  is  the  minimum  capacity  cut. 

For  network  1  all  minimal  cuts  Y  and  conditional  performance 
measures  are  tabulated  in  Table  VIII.  As  the  number  of  minimal 
cuts  in  network  2,  3  and  4  is  large  (17,  15  and  78  respectively) 

they  are  first  ordered  by  decreasing  criticality  indices  and  only 
first  ten  are  tabulated  in  Tables  IX,  X  and  XI  respectively. 
Algorithm  by  Tsukiyama  et  al  [1980]  is  used  to  enumerate  all 
minimal  (s,t)  cuts  in  a  network. 

The  algorithms  proved  to  be  very  efficient  and 
computationally  stable.  It  took  about  one  second  of  CPU  time  to 
compute  the  distribution  of  the  maximum  flow  for  each  network  and 
less  than  0.2  seconds  per  cut  to  evaluate  the  conditional 
performance  on  IBM  4341  system.  Using  simulation  to  obtain 
equally  accurate  estimates  would  no  doubt  take  considerably 


longer  time. 


8.  Conclusion. 


We  have  provided  a  method  of  computing  the  exact 
distribution  of  the  maximum  flow  in  planar  networks  with 
exponentially  distributed  arc  capacities  by  constructing  a 
continuous  time  Markov  chain  with  one  absorbing  state  so  that  the 
time  until  absorption  in  this  state  is  equal  to  the  maximum  flow 
from  source  to  sink  in  the  network.  The  state  space  of  this  CTMC 
is  shown  to  be  the  set  of  all  paths  in  the  network,  augmented  by 
an  absorbing  state.  When  the  paths  are  put  in  "topmost  path 
first”  order  the  rate  matrix  of  this  chain  becomes  upper 
triangular.  This  makes  all  calculations  very  straightforward. 
Algorithms  are  developed  to  compute  the  distribution  and  moments 
or  the  maximum  flow. 

Further,  a  modified  CTMC  is  constructed  to  compute  the  joint 
distribution  of  the  criticality  indicator  of  a  given  minimal 
( 3 , t  )  cut  and  the  maximum  flow.  Algorithms  are  developed  to 
compute  the  probability  that  the  given  cut  is  a  minimum  capacity 
cut  in  the  network. 

Note  that  it  is  possible  to  convert  a  maximum  flow  problem 
for  a  planar  network  with  exponential  arc  capacities  into  a 
shortest  path  problem  in  the  dual  network  with  exponential  arc 
lengths  and  solve  this  shortest  path  problem  by  the  method 
described  in  Kulkarni  [1984],  There  are  two  disadvantages  to 
this  approach.  First,  one  has  to  construct  the  dual, 
considerable  extra  work  for  large  networks.  Second,  the  state 
space  requirement  of  this  equivalent  shortest  path  problem  is 
almost  always  more  than  that  of  the  path  filling  algorithm 


described  here.  It  is  important  to  solve  the  maximum  flow 
problem  directly  since  it  provides  many  interesting  algorithms 
and  concepts. 

One  major  drawback  of  CTMC  method  presented  here  is  that  the 
size  of  the  state  space,  which  is  equal  to  the  number  of  paths, 
can  grow  exponentially  with  network  size.  Thus  the  time  and 
space  requirements  of  the  algorithms  can  increase  rapidly.  It 
should  be  noted,  however,  that  for  very  large  networks,  any  other 
method  of  analysis  will  also  be  slow. 

An  extension  of  our  method  to  non-planar  networks  does  not 
seem  possible,  since  the  "topmost  path  first"  order  does  not  make 
sense  in  non“planar  networks.  It  seems  to  us  that  an  entirely 
different  approach  would  be  needed  to  solve  the  maximum  flow 
problem  in  non-iplanar  networks. 

A  quantity  that  io  of  interest  in  network  flow  analysis  is 
the  utilization  of  an  arc,  i.e.  the  fraction  of  the  capacity 
utilized  when  the  flow  in  the  network  is  maximum.  This  and  some 
related  quantities  will  be  studied  in  a  forthcoming  paper. 
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Table  III,  The  Y-admissible  Paths 
in  Network  1.  Y  «  {3,5,9}  - 


i 

Nodes 

visited 

by  path  i 

1 

1 

2 

4 

5 

2 

1 

2 

3 

4  5 

3 

1 

2 

3 

6  5 

4 

1 

3 

4 

5 

5 

1 

3 

6 

5 

Table  IV.  The  Q(Y)  Matrix  for  Network  1 
X  -  13,5,9} 

1 _ 2 _ 3 _ 2 _ 5 _ 6  7 


u(Pl  ) 

0 

0 

0 

0 

0 

-u<P2) 

U5 

Ml  +U4 

0 

0 

0 

0 

- u ( P3  ) 

0 

Ml  +M4 

M9 

0 

0 

0 

" M ( P4  > 

^5 

0 

0 

0 

0 

0 

~  M ( P5  ) 

M9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

State  6  is  the  "good"  absorbing  state. 
State  7  is  the  "bad"  absorbing  state. 


I 


uj  r\) 


X 

Network  1 

Network  2 

Network  3 

Network  4 

-3 

0.0 

0.0 

0.0 

0.0 

-2 

0.0 

0 . 0 

0.0 

0.0 

-  1 

0.1391 

0.1421 

0.1 354 

0.1487 

0 

0.5727 

0.5680 

0.5788 

0.5568 

1 

0.8502 

0.8496 

0.8506 

0.8473 

2 

0. 9566 

0.9574 

0.9554 

0.9596 

3 

0 . 9887 

0.9893 

0 . 9880 

0.9910 

Table  VIII.  Conditional  Performance:  Network  1. 


Cut  # 

Y  -  C(X,X) 

R(Y) 

u  (Y) 

o  ( Y  ) 

1 

3.5.6 

0.2678  “ 

T7oTo9 

0.5787 

2 

1,2 

0. 1780 

0.8949 

0.5632 

3 

2.3.4 

0. 1 608 

1 .0461 

0.5765 

4 

6.7.8 

0.1382 

1 .0554 

0.5785 

5 

7.9 

0.1142 

0.9056 

0.5657 

6 

1.5,6 

0.0680 

0. 9989 

0.5683 

7 

3,5.9 

0.0535 

1  .0056 

0.5698 

8 

1.5.9 

0.0136 

0.9636 

0.5588 

9 

2. 4, 7. 8 

0.0059 

1.0138 

0.5154 

R(Y) 

u(Y) 

! 

o(Y) 

0.2204 

0.9776 

0.6049 

0.1405 

1  .0833 

0.6060 

0.0923 

1.1483 

0.6263 

0.0878 

1  .0833 

0.6060 

0.0798 

1  . 3066 

0.6409 

0.0770 

1  .3066 

0.6409 

0.0722 

1.2394 

0.6164 

0.0499 

1  . 3066 

0. 6409 

0.0481 

1  .  2990 

0.6393 

0.0451 

1 .2394 

0.6164 

Table  X.  Conditional  Performance:  Network  3. 


Cut  » 

Y  -  C(X,X) 

R(Y) 

u  ( Y ) 

o  ( Y  ) 

1  .2 

0.3051 

r  0.9899  ' 

0.6516 

1.6 

0. 1566 

0.8467 

0.5590 

7,8,11 

0.1045 

1 .0708 

0.6292 

3,8,11 

0.0871 

1 .0708 

0.6292 

1.11 

0.0687 

0.7289 

0.4947 

4,6,7 

0.0621 

1.1128 

0.6337 

6, 7, 8, 9 

0.0529 

1.3103 

0.6712 

8 

3, 6. 8, 9 

0.0441 

1.3103 

0.6712 

mam 

4,6,7 

0.0211 

0. 9790 

0.5807 

3.4,11 

0.0175 

0 . 9790 

0.5807 

i 

Table  XI.  Conditional  Performance:  Network  4. 


Figure  1.  The  Example  Network  1. 

The  numbers  on  the  arcs  are  arc  numbers. 

The  numbers  in  brackets  represent  mean  capacities 
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